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Exact analytical solutions of the time-dependent Schrodinger equation with the initial condition 
of an incident cutoff wave are used to investigate the traversal time for tunneling. The probability 
density starts from a vanishing value along the tunneling and transmitted regions of the potential. 
At the barrier width it exhibits, at early times, a distribution of traversal times that typically has 
a peak rp and a width Ar. Numerical results for other tunneling times, as the phase-delay time, 
fall within Ar. The Biittiker traversal time is the closest to Tp. Our results resemble calculations 
based on Feynman paths if its noisy behaviour is ignored. 

PACS numbers: 03.65.Bz, 03.65. Ca, 73.40. Gk 
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Quantum tunneling, that refers to the possibility that 
a particle traverses through a classicaUy forbidden region, 
constitutes one of the paradigms of quantum mechanics. 
In the energy domain, where one solves the stationary 
Schrodinger equation at a fixed energy E, tunnehng is 
well understood. In the time domain, however, there are 
aspects still open to investigation. Recent technological 
achievements as the possibility of constructing artificial 
quantum structures at nanometric scales 0] or the ma- 
nipulation of individual atoms have stimulated work 
on time-dependent tunneling both at applied and funda- 
mental levels. A problem that has remained controversial 
and the subject of a great deal of attention over the years 
is the tunneling time problem ||^, that can be stated as 
the question: How long it takes to a particle to traverse a 
classically forbidden region? In time-dependent tunnel- 
ing, many works attempting to answer the above question 
consider the numerical analysis of the time-dependent 
Schrodinger equation with the initial condition of a Gaus- 
sian wave packet . A common feature of the ma- 
jority of these approaches is that the initial wave packet 
extends through all space. As a consequence the initial 
state, although it is manipulated to reduce as much as 
possible its value along the tunneling and transmitted 
regions, contaminates from the beginning the tunneling 
process and hence usually it is required a long time anal- 
ysis of the solutions. The above situation may be cir- 
cumvented by considering cutoff wave initial conditions 
0,||jl||l|,|l§. 

In this work we consider a analytic time-dependent so- 
lutions to the Schrodinger equation with the initial condi- 
tion at r = of a incident cutoff wave, to investigate the 
traversal time for tunneling through a potential barrier. 



The problem may be visualized as a gedanken experiment 
consisting of a shutter, situated at a; = 0, that separates 
a beam of particles from a potential barrier of height Vq 
located in the region < x < L. AtT = the shutter 
is opened. The probability density rises initially from a 
vanishing value and evolves with time through x > 0. At 
the barrier edge x — L, the probability density at time r, 
yields the probability of finding the particle after a time 
T has elapsed. Since initially there is no particle along 
the tunneling region, detecting the particle at the barrier 
edge at time r, provides a measure of its traversal time 
through the tunneling region. 

The transient behavior of the time-dependent solution 
at early times and at distances close to the interaction 
region plays a significant role in our approach. Other 
formulations, based on the stationary solutions of the 
Schrodinger equation [^j^, refer to asymptotically long 
times at large distances and hence ignore transient ef- 
fects. These approaches provide a single value for the 
traversal time. In contrast, our approach leads to a distri- 
bution of traversal times as in works based on the Feyn- 
man path integral method though as indicated 
below both approaches differ in important aspects. 

In a recent paper we have obtained the time-dependent 
solution to the Schrodinger equation for tunneling 
through an arbitrary potential of finite range with the 
initial condition of a cutoff plane wave of momentum k. 
The solution may be written as a term proportional to 
the free solution plus a contribution involving an infinite 
sum of resonance terms associated with the S*— matrix 
poles of the potential . Our approach is based on the 
Laplace transform technique and considers some analyt- 
ical properties of the outgoing wave propagator. Some 
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decades ago Moshinsky considered the free case solution 
to the above problem Moshinsky showed that the 

probability density, for a fixed value of the distance xq 
as a function of t, exhibits a transient regime that he 
named diffraction in time. Recently, observations of that 
phenomenon have been reported For the sake of 

simplicity in our approach, as Moshinsky also did, we 
consider the instantaneous removal of the shutter. This 
may be seen as a kind of 'sudden approximation' to a 
shutter opening with finite velocity, where the treatment 
becomes more involved [Q. As shown below the terms 
depending on the S*— matrix poles provide a novel tran- 
sient behavior that may dominate the early times in the 
tunneling process. 

The plane wave cutoff initial condition discussed in 
Refs. ||ll|,|l^ refers to a shutter that acts as a perfect 
absorber (no reflected wave). One can also envisage a 
shutter that acts as a perfect reflector. In such a case 
the initial wave may be written as, 



ipsix,k,T = 0) 



-ikx 



0, 



X < 
a; > 0. 
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One can then proceed along lines similar to those dis- 
cussed in Ref. to derive the time-dependent solution 
ips{x,k,T) of the Schrodinger equation for a potential 
V{x) that vanishes outside the region < x < L. The 
solution along the internal region reads, 

ipsix, k, t) = (j){x, k)M{0, k, t) - (t>{x, ~k)M{0, -k, r) 

OO 

-^0„(x)M(O,A:„,t), (0<x<i) (2) 

n 

where (j){x, k) refers to the stationary solution and 
4>n{x) — 2ikun{0)un{x) / {k'^ — fc^) is given in terms of 
the resonant (Gamow) states {w„(x)} and complex poles 
{kn} of the problem |16 1^. Similarly the transmitted 
solution jl5| becomes, 

ijs{x,k,T) ^T{k)M{x,k,T) -T{-k)M{x,~k,T) 

OO 

-^T„M(a:,fc„,r), {x > L) (3) 

n 

where T{k) and T{—k) are transmission amplitudes, and 
r„ = 2ikun{Q)un{L)ex-p{—iknL)/{k'^ the above 

two equations the functions M(x,k,T) and M(a:,fc„,T) 
are defined as, 



Mix,q,t) ^ ^e^'^-"/^''^^ey"erk{yg), 



where the argument yq is given by 
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In Eqs. (^) and (|^) q stands either for k oi kn, the index 
n refers to a given complex pole. Poles are located on the 



third and fourth quadrants of the complex fc-plane. The 
solution for the free case with the reflecting initial condi- 
tion is, V's (2;, fc, t) = M{x, fc, t)—M (x, — fc, t). ^From the 
analysis given in Ref. ||l^ one can see that the above ex- 
act solutions satisfy the corresponding initial conditions, 
i.e., they vanish exactly for x > 0. At very long times it is 
also shown in Ref. that the terms M{x, kn, r) that ap- 
pear in the above equations vanish. The same occurs for 
M{x, —k,T) while, as shown firstly in Ref. , M{x,k,T) 
tends to the stationary solution. Hence, at long times, 
each of the above exact solutions go into the correspond- 
ing stationary solutions, namely, along the internal region 
as ^j{x,t) = (f>(x,k)exp{—iET/h) and along the exter- 
nal region as ^/'(a;, r) = T(k)e'xp{ikx)e'xp{—iET/h). Note 
that at early times and short distances there is a com- 
petition between the contribution of the free-type terms 
{M functions depending on k) and the pole terms {M 
functions depending on either fc„ or fc_„) in Eqs. (H) and 
(^. As exemplified below, depending on the potential 
parameters one may have the predominance of one or 
the other type of terms. Note also that the initial state 
is not strictly monochromatic (it extends from —00 to 0) 
and hence it has a distribution of components around k 
in momentum space. One could construct an initial cut- 
off wavepacket as a linear combination of cutoff waves. 
However, since we compare below with definitions of tun- 
neling times involving plane waves, wavepackets will not 
be considered here. Besides, in general they involve no 
negligible momentum components above the barrier po- 
tential and hence obscure the dynamics of tunneling. 

In order to apply the above ideas, we consider a model 
that has been used extensively for the tunneling time 
problem, namely, the rectangular barrier potential, char- 
acterized by a height Vq in the region < a; < L. The 
shutter is located at a; = 0. In order to calculate Eqs. (|^) 
and (III) for the initial condition (0), in addition to the pa- 
rameters Vq, L, and that corresponding to the incident 
energy E = h^k^ /2m, we need to determine the complex 
poles {kn} and resonant states {u„(x)}. It is well known 
that for a finite range potential there are an infinite num- 
ber of poles. The 5— matrix poles for the rectangular 
barrier potential may be obtained from the correspond- 
ing transmission amplitude T{k) = Akqexp{—ikL)/ J{k), 
where q = [k^ — kg]^/^ with /cq = 2mVo/h^. They corre- 
spond to the zeros of J{k) in the fc— plane, namely, 

J{k) = {q + k)^cxp{-iqL) - {q - k)^exp{iqL) = 0. (6) 

We follow a well established method to obtain the solu- 
tions to the above equation jl^,^. The resonant states 
of the problem satisfy the time-independent Schrodinger 
equation of the problem with outgoing boundary condi- 
tions They read, 

wn(x) = C„ [e'«"" + 6„e-^'""] , (0 < a; < L) (7) 

where &„ = (g„ -f kn)/ {qn — kn) and C„ may be obtained 
from the normalization condition p^ . 
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Note that both the complex poles and the resonant states 
are a function of Vq and L and hence are a property of 
the system. 

To exemplify the time evolution of the probability den- 
sity we consider the set of parameters: Vb = 0.711 e^, 
L = 10 nm, E = 0.1422 eV, m* = 0.067 We, inspired in 
semiconductor quantum structures ^ . Our choice of pa- 
rameters guarantees that most momentum components 
of the initial state tunnel through the potential. Fig- 
ure |l| shows a plot of |-0(L, r)p, calculated at the barrier 
edge X = L as 'A function of time in units of the free pas- 
sage time Tf = mL/{fik) = 11.56 /s. We have used Eq. 

though Eq. @ holds the same. The time-dependent 
solution is normalized by |T(fc)p = 5.332 x 10~^. One 
sees that as soon as r 7^ the probability density starts 
to grow up. As discussed elsewhere this is due to 
the non-relativistic character of the description. Einstein 
causality may be fulfilled by cutting off the contributions 
to the probability density smaller than tq = L/c. In 
our example ro = 0.033 /s or to/t/ = 0.0028, too small 
to be appreciated in Fig. |l|. At early times one sees a 
time domain resonance structure. Thereafter the proba- 
bility density approaches essentially its asymptotic value. 
We found that the resonant sum is the relevant contri- 
bution to the time domain resonance since that of the 
free-type term is quite small and varies smoothly with 
time. In the transmitted region, x > L, not shown 
here, the time domain resonance becomes a propagat- 
ing structure, as follows from Eq. (|^). The time do- 
main resonance corresponds to a transient effect and as 
it propagates through the transmitted region becomes 
smaller and smaller. Asymptotically, at large distances 
and times, it becomes very small while the free-type term 
becomes the dominant contribution with its wavefront 
propagating with velocity v = hk/m. Calculations using 
the absorbing initial condition exhibit a similar time do- 
main resonance. Hence a linear combination of reflecting 
and absorbing initial conditions should also exhibit it. 

The time domain resonance represents a distribution of 
traversal times. The corresponding peak represents the 
largest probability to find the tunneling particle at the 
barrier edge x — L. In our example, as shown in Fig. 
the time domain resonance peaks at Tp = 5.326 /s, faster 
than the free passage time across the same distance of 
10 nm, that is, Tp/r/ = 0.46. Note that the distribu- 
tion is quite asymmetric. Altough the first resonance 
term of the solution provides the main contribution, con- 
vergence of the series usually requires to sum up to 100 
terms. The inset displays the probability density from 
t/t/ — 2 up to r/r/ — 20. One sees a small struc- 
ture around t/t/ = 3 and then the probability density 
decreases very fast towards unity, the stationary regime. 
The main range of traversal times occurs around the peak 



value Tp. We define the width of the distribution, Ar, by 
the rule of the half-width at half-maximium. This yields 
Ar = 13.48 /s or Ar/r/ = 1.16. The resonance is broad, 
since At « 2Tp. We have found that for fixed Vb and 
E, and a decreasing L, the width diminishes. The same 
occurs for fixed E and L, and an increasing Vq. System- 
atically, however, Ar > Tp. For the sake of comparison, 
the arrows in Fig. ^ indicate the values calculated for a 
number of definitions of tunneling times existing in the 
literature for the rectangular barrier potential |22j, as 
the Larmor time of Baz and Rybachenko, tlm', the semi- 
classical or Biittiker-Landauer time, tbl', the Biittiker 
traversal time, rs, and the phase-delay time, tjj [^ . 
All of them fall within the broad range of values given by 
Ar. Note, however, that the Biittiker traversal time tb 
is the closest to Tp. As shown below we have found this 
situation extensively in our numerical calculations. Also, 
since the barrier is opaque, tb is close to tbl- 

We refer briefly to approaches to the tunneling time 
problem based on the Feynman path integral method 
For plane waves and a rectangular barrier po- 
tential Fertig has derived an expression, C(r), that 
gives the probability amplitude that a particle remains 
a time r in a region, (Eq. (3) of Ref. 0). Recently 
Yamada [0 has plotted G{t) = |C(r)l^ versus r (Fig. 
2 of Ref. 13). His parameters are the same as in our 
Fig. 1^, i.e., Vo/E = 5 and kL = 5. Our calculation 
for |?/'(L,t)P resembles the average shape of G(r), pro- 
vided its noisy behavior is ignored. Note, however, that 
the meaning of both quantities is different. As indicated 
by Yamada, G(r) refers to a 'residence time' j2^] whereas 
our approach corresponds to a 'passage' or traversal time 

In Fig. 1^ we plot Tp (solid squares)for different values 
of the opacity a = koL, with fcp = [2mVb]^/^/?i. Keeping 
Vq fixed and varying L defines a{L). We can then identify 
two regimes, one in the range 2 < a{L) < 5, the tunnel- 
ing regime, where Tp remains almost constant as a{L) 
increases, and another regime, the opaque regime, with 
a{L) > 5, where we find that Tp increases linearly. The 
first behaviour above is related to the first top-barrier 
S-matrix pole and the second one to the components of 
the incident wave that go above the barrier. There is still 
another regime, not shown in Fig. ^, where a{L) < 1, 
that corresponds to very shallow or very thin barriers or 
both and will not be considered here. There the free-type 
terms in Eq. (^ dominate over the resonant contribu- 
tion. For comparison we plot the Biittiker traversal time 
Tb (hollow circles). We see that tb remains rather close 
to Tp. Note, however, that tb behaves linearly in the 
whole range. This different qualitative behaviour as a 
function of L between both times deserves further study. 
The inset in Fig. ^ exhibits a similar comparison for the 
opacity a(Vo), with L fixed and varying Vq. Here we 
observe that tb remains quite close to Tp in the whole 
range. Regarding the phase-delay time td, its predic- 



3 



tions usually fall within the width At. For fixed Vq and 
i?, td as a function of L exhibits qualitatively a different 
behaviour than that of Fig. |l| (See Fig. 5 in ref. ||]). 

To end we stress that the largest probability to find 
the tunneling particle at the barrier width, given by r^, 
is sensitive to both variations of the barrier width L and 
of the height Vq, and also, that the Biittiker traversal 
time is found very close to the value of Tp though we find 
qualitative differences between them as a function of L. 
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FIG. 1. Plot of |i/)(L,t)|^ at the barrier edge 
X = L = IQ nm as a function of time in units of the free 
passage time r/. The inset shows \ip{L,T)\'^ at larger times. 
The arrows indicate the values of the Larmor time LM, the 
semi-classical time BL, the Biittiker traversal time B, and the 
phase-delay time D. See text. 

FIG. 2. Plot of the exact time domain resonance peak Tp 
(solid squares) versus the opacity a(L) (Vb fixed). For com- 
parison we plot the Biittiker traversal time tb (hollow circles). 
The inset shows a similar calculation versus the opacity a(Vb) 
(L fixed). See text. 
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